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SECTION  I 


INTRODUCTION 

For  quite  some  time,  there  has  been  interest  in  the  penetration  of 
semi-infinite  targets  by  long  rods  and  jets.  Numerous  experimental  and 
numerical  studies  of  the  problem  have  been  carried  out  over  the  past  four 
decades.  The  numerical  studies  have  now  evolved  into  gigantic  hydrocodes 
which  attempt  to  account  for  every  detail  of  the  event.  These  codes  are 
expensive  to  run  and  do  not  offer  much  engineering  insight  into  the  pene¬ 
tration  process.  It  is  for  this  reason  that  the  one-dimensional,  eroding- 
rod  model  of  Tate  (References  1  and  2)  is  both  a  popular  and  useful  tool  for 
describing  the  event.  For  nearly  two  decades,  Tate's  theory  has  been 
regarded  as  the  foundation  for  simple  engineering  modeling  of  the  penetra¬ 
tion  process. 

In  this  paper,  a  generalization  of  the  basic  equation  of  motion  of 
the  rod  is  made  which  leads  to  a  modification  of  Tate'3  theory.  The 
modified  theory  accounts  for  the  expansion  of  the  penetrator  tip  into  a 
mushroom  of  diameter  greater  than  that  of  the  original  rod,  and  provides  an 
accurate  treatment  of  mass  transfer  from  the  undeformed  segment.  Compari¬ 
son  with  Tate's  theory  is  made  when  appropriate. 


SECTION  II 


THEORY 

1 .  GENERAL 

The  approximate,  one-dimensional  analysis  of  the  penetration  of  a 
half-space  by  a  long  (slender)  rod  is  best  explained  for  the  case  in  which 
the  penetration  is  accompanied  by  consumption  of  the  rod.  The  other  two 
cases,  impact  without  penetration  and  rigid-body  penetration,  can  be  con¬ 
sidered  as  special  cases  of  the  more  general  situation  of  the  penetrator 
being  consumed  as  it  penetrates. 

Figure  1  shows  the  penetrator,  a  cylindrical  rod  of  initial  length  L, 
density  p,  and  cross-sectional  area  A.  Let  V  denote  the  initial  velocity 
with  which  it  impacts  the  semi-infinite  target.  After  some  time  t  has 
elapsed,  the  rod  has  penetrated  the  target  to  some  depth  z.  In  the  process, 
a  portion  of  the  rod  of  length  X  has  been  consumed.  The  remaining  portion 
of  the  rod  has  length  L-X,  denoted  by  £;  it  is  assumed  to  move  as  a  rigid 
body  with  remaining  velocity  v.  The  speed  of  penetration  into  the  target, 
z,  is  denoted  by  u.  Here  the  superposed  dot  indicates  derivative  with 
respect  to  time. 

Figure  2(a)  shews  the  rigid  end  of  the  rod  at  time  t.  An  instant 
later,  at  time  t  +  At,  a  portion  of  the  rod  end  has  been  consumed  by  the 
penetration  process  and  the  remaining  rod  end  has  been  decelerated.  This 
is  depicted  in  Figure  2(b).  From  these  figures  a  simple,  one-dimensional 
impulse-momentum  equation  can  be  written.  Because  the  internal  forces  F  are 
equal  and  opposite,  only  the  external  force  P  contributes  to  the  impulse. 

As  it  is  opposite  to  the  directio.,3  of  u  and  v,  its  contribution  is  negative 
-PAt.  The  total  momentum  at  t  ♦  At  is  pAAXu  +  pA(L-X~AX)(  v  +  Av).  At  t  the 


Figure  1.  Schematic  of  Hod:  (a)  Shows  Plastic  Portion  X,  and 
Undeformed  Portion  L  ~  X  =  l;  (b)  Shows  Penetration 
into  Target  to  a  Degth  z. 


momentum  was  pA(L-X)v.  Subtracting  to  obtain  the  momentum  change  and 
equating  that  to  the  impulse  gives 

-PAt  -  pA{(L-X)Av-AX(v-u)-AXAv}  (1) 

Dividing  through  by  At  and  taking  limits  as  At  approaches  zero  gives 

-P  *  pAKL-X)v-X(v-u)}  (la) 

as  the  equation  of  motion  for  the  rigid  end  of  the  rod.  From  the 
definition  d  -  L-X,  Equation  (la)  can  be  recast  in  the  form 

-P  =  pA  (d’v+"d(  v-u)}  Ob) 

Again  referring  to  Figure  2(b),  it  is  assumed  that  a  reasonable  approxima¬ 
tion  to  the  external  force  is  given  by 

P  =  pA/(1+e)  (2) 

Here  e  denotes  the  engineering  strain  in  the  deformed  element  of  the  rod 

so  that  A/(1+e),  assuming  isochoric  deformation,  is  its  current  cross- 
sectional  area;  and  p  denotes  the  interface  pressure  at  the  penetrator  tip. 
Substitution  of  Equation  (2)  into  Equation  (1b)  gives 

dv  +  d(v-u)  -  -p/p(1+e)  (3) 

This  is  the  rod  equation  of  motion. 

With  reference  to  Figure  2,  during  the  time  interval  At  the  back  end 
of  the  remaining  penetrator  moves  a  distance  vAt  while  its  front  end  moves 
uAt .  The  consequent  change  of  length  Ad  is  (u-v)At  so  that 

d  *  -(v-u)  (J») 

This  is  a  purely  kinematical  relationship  among  the  variables. 

To  complete  the  statement  of  the  problem  requires  the  specification  of 
functional  relations  for  u,  p,  and  e.  Suppose,  for  example,  that  the 
penetration  velocity  u  is  some  function  of  the  rod  speed  v,  say  u-f(v).  It 
is  frequently  assumed  that  the  pressure  p  is  a  function  of  penetration 
velocity;  in  this  case  p  is  also  a  function  of  v:  p  =  p(u)  =  p(f(v))  =*  g(v). 
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The  strain  e  can  be  assumed  to  depend  upon  pressure  (or  stress)  p;  then  it 
too  is  a  function  of  v:  e  -  e(p)  «  e(g(v))  -  h(v).  The  rod  equation  of 
motion  can  now  be  written  as 

8,v  -  [v-f(v)]2  -  ~g(  v)/{p[1 +h(v)  ])  (5) 

Equation  (5)  can  be  solved  generally  because  the  variables  are  separable. 
By  making  the  substitution  v  »  -(dv/d£)[  v-f( v)]  and  rearranging,  Equation 
(5)  becomes 


-£(dv/dJl)  -  [v-f  ( v)  ]~g(  v) / {p[1  +h(  v)  ][v-f  ( v)  ]} 
Separating  variables  and  integrating  from  the  starting  conditions  i 
v  «  V  leads  to 


inU/L) 


dv 


g(v)  {p[l+h(v)][v-f(v)]}_1  -  [v-rf(v)] 


(6) 

-  L  and 


(7) 


This  is  the  general  dependence  of  remaining  rod  length  upon  remaining  rod 
speed. 

Several  special  cases  can  be  derived  from  the  foregoing  general 
theory.  For  example,  the  Tate  analysis  (References  1  and  2)  is  obtained  by 
assuming  that  (v-u)  and  e  are  negligible  in  Equation  (3),  and  that  p  is 
given  by  a  modified  Bernouli  equation 

p  -  (1/2)p2pu2  +  R  -  (1/2)p(v-u)2  ♦  Y  (8) 

Here  p2p  denotes  target  density,  and  R  and  Y  are  the  respective  strength 

factors  of  target  and  penetrator.  In  general  these  can  be  quite  different 
from  the  static  yield  strengths  of  the  two  materials.  Strain  rate  effects 

tend  to  increase  these  strengths  while  thermal  softening  tends  to  decrease 

them.  Equation  (8)  is  subject  to  the  restrictions  that:  the  right  hand 
equality  is  invalid  when  the  rod  penetrates  as  a  rigid  body;  both  are 
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invalid  when  there  is  no  penetration.  Unsteady  flow  is  entirely  neglected 
in  Equation  (8)  and  at  present  there  is  no  satisfactory  way  of  estimating 
the  significance  of  this  omission. 

In  this  paper  a  more  general  form  of  the  Tate  analysis  will  be 
treated.  It  is  assumed  that  Equation  (8)  gives  the  stagnation  point 
pressure  with  sufficient  accuracy,  but  (v-u)  and  e  are  retained  in  Equation 
(3).  Inasmuch  as  Equation  (8)  arises  from  pseudo-steady-state  considera¬ 
tions,  it  seems  commensurate  to  assume  that  e  is  simply  constant.  This 
assumption  is  made  for  mathematical  convenience  and  from  ignorance.  While 
it  is  expected  that  the  appropriate  value  for  e  will  grow  larger  in  some 
manner  with  V,  the  details  of  any  such  relationship  remain  to  be  developed. 
However,  treating  e  as  some  constant  in  any  given  impact  is  in  keeping  with 
the  other  steady-state  assumptions. 

Combining  Equations  (3),  (4),  and  (8)  gives  the  rod  deceleration  as 
either 

v  »  [(1  +e)  (v~u)2  -  (1/2)u2u2  -  A2c2]/Sl(1  +e)  (9) 

or 

v  *  [(1 /2+e) (v-u)2  -  c2]/i(1+e)  (9a) 

Here  A2  denotes  the  ratio  R/Y  and  c2  the  quantity  Y/p  which  has  the 
dimensions  of  velocity  squared.  Equations  (*0,  (8),  and  (9)  or  (9a) 
comprise  a  system  of  three  equations  in  three  unknowns:  u,  v,  and  2,.  In 
the  process  of  solving  them,  however,  another  equation 

z  «  u  (10) 

should  be  utilized  in  order  to  find  the  most  important  element  in  the 
problem,  viz.,  the. penetration  depth.  The  problem  now  is  completely 
formulated.  The  parameters  embedded  in  it  are  conveniently  taken  as  u2  » 
(target  density /penetrator  density),  A2  »  (target  strength/penetrator 


strength),  c2  -  Y/p,  e  »  engineering  strain  in  the  deformation  zone  of  the 
penetrator,  and  V  -  the  initial  speed  of  impact  of  the  rod.  No  general 
solution  beyond  Equation  (7)  is  available  for  the  foregoing  set  of  equa¬ 
tions;  they  must  be  integrated  numerically  for  particular  values  of  the 
parameters. 

Before  any  calculations  are  attempted,  the  possibilities  mentioned 
earlier,  failure  to  penetrate  and  penetration  without  rod  consumption,  must 
be  considered.  Intuitively,  the  former  would  be  associated  with  large 
values  of  p2  ,  large  values  of  X2  and/or  small  values  of  V;  and  the  latter 
would  correspond  to  p2  and/or  X2  being  small. 

2.  NO  PENETRATION 

Equation  (8)  can  be  used  to  find  u  in  terms  of  v,  u  ,  and  c  . 

Negative  values  of  u  imply  no  penetration.  In  this  case  set  u  -  0  and  p  * 
(1/2)pv2  +  Y.  Then  from  Equation  (9a) 

v  *  [ ( 1 +e) v2-x2c2]/£( 1 +e)  (11) 

From  Equation  (U)  v  »  -£  and  v  -  -j,” which  can  be  substituted  into  Equation 
(11)  to  obtain 

IZ'+  i2  -  X2c2/( 1 +e)  (12) 

The  solution  of  Equation  (12),  taking  account  of  the  initial  conditions 
l  =  L  and  i  ■  -V  at  t  •  0,  is 

lr/L  -  [1 ~V2(1 +e)/X2c2]1 /2  (13) 

where  if  denotes  the  length  of  rod  remaining  after  it  comes  to  rest.  For 
initial  velocities  large  enough  that  V2(1+e)  >  X2c2  Equation  (13)  will  not 
apply.  Instead  the  rod  will  be  entirely  consumed  against  the  face  of  the 


target. 


3.  NO  ROD  DEFORMATION 


Values  of  u  larger  than  v  imply  rigid  body  penetration.  In  this  case 
set  u  -  v  and  p  -  (1/2)  y2pu2  +  Then  from  Equation  (9) 

A  »  -[X2c2+(1/2)y2u2]/i(1+e)  (IK) 

Now  Equation  (4)  shows  that  i  -  0  so  that  l  remains  constant.  Noting  that 
<5  *  u(du/dz)  Equation  (14)  can  be  readily  integrated  to  obtain 

z  «  zQ  +  [!0( 1 +e) /u2]£n  [(2X2c2  +  y2u2)/(2X2c2  +  y2u2)]  (15) 

Here,  following  Tate  (Reference  2),  u0  denotes  the  rod  velocity  at  which  u  - 
v  first  occurs  and  %0  is  the  corresponding  length.  Also,  z0  denotes  the 
penetration  to  this  point.  The  increment  of  penetration,  ze,  that  occurs  in 
bringing  the  rigid  rod  to  rest  is  obtained  by  substituting  u-0  into  Equation 
(15).  Thus,  the  total  penetration  depth  is  given  by 

Zf  -  zQ  +  [l0(1+e)/u2]lln[1  ♦  y2u2/2x2c2]  (16) 

Note  that  this  result  is  almost  identical  to  Tate's  and  has  the  general  form 
of  the  Petry  equation  (Reference  3)  or  the  DeMarre  formula  (Reference  4). 

For  the  situation  in  which  rigid  body  penetration  proceeds  from  initial 
impact,  uQ  -  V,  !l0  »  L,  and  Equation  (16)  yields 

zf  -  [L(1  +e)/y2]!ln[1  +y2V2/2X2c2]  (17) 

4.  PENETRATION  ACCOMPANIED  BY  ROD  DEFORMATION 

Having  disposed  of  the  rigid-body  and  zero  penetration  cases,  the 
penetration-with-consumption  case  is  resumed.  The  value  of  u  obtained  from 
Equation  (8)  depends  upon  the  parameters  y2  and  X2.  When  either  of  these 
is  unity,  the  equation  is  essentially  linear;  otherwise  it  is  quadratic  and 
rather  more  complex.  Consider  first  the  condition  of  equal  strengths  of 
target  and  rod;  here  R  =  Y  so  that  x2  -  1  and  Equation  (8)  yields 


In  this  situation  u  is  always  smaller  than  v  and  the  two  reach  zero 


simultaneously.  There  is  no  rigid-body  penetration  phase  under  these 
conditions.  From  Equation  (18)  (v-u)  -  yv/(1+p)  so  that  Equation  (H) 
can  be  written  as 

l  -  -pv/(1+ji)  (19) 

Differentiating  Equation  (19)  and  substituting  into  Equation  (9a)  gives 

li'+ai2  -  b  -  0  (20) 


where 

a  =*  u(  1  /2+e)/(  1  +p) ( 1  +e)  (21) 

and 

b  -  uc2/(1+p)(1+e)  (22) 

By  defining  5  =  £2,  from  which  d£/d£  »  2St’,  Equation  (20)  can  be  rewritten 
as 

£(d£/d£)  +  2a£  -  2b  -  0  '  (23) 

/ 

This  equation  is  linear  and  can  be  simply  integrated  subject  to  £(L)  • 
u2V2/(1*u)2. 

Two  cases  develop:  first,  a-0,  which  implies  e--1/2.  In  this  case 
Equation  (23)  integrates  to 

E,  «  £2  -  ( y2V2) / ( 1  +u )2  -  2b£n(L/£)  (2*4) 

From  Equation  (2M)  it  can  be  seen  that  there  will  always  exist  a  non-zero, 
unconsumed  rod  length  If  corresponding  to  £  -  0  and  given  by 

£f  -  L  exp{-y2V2/2b(1 ♦u)2}  -  L  exp{-yV2/c2(1 +p)}  (25) 

For  this  case  in  which  u  is  directly  proportional  to  v,  Equation  (18), 
i  is  also  proportional  to  v,  Equation  (19).  From  Equations  (18)  and  (19) 
u  »  -i/u  so  that  the  penetration  depth,  given  by  Equation  (10),  in  this  case 
becomes 

ft  ft 


The  total  depth  of  penetration,  zf,  is  given  by 


zt  -  (L-lf)/u  (26a) 

where  If  is  given  in  Equation  (251. 

The  more  important  case  is  a  *  0.  For  this  case,  it23-1  is  an 
integrating  factor  of  Equation  (20)  and  the  solution  is 

6  -  i2  -  b/a  +  [y2V2/(1+y)2  -  b/a] (L/l)2a  (27) 

If  a  is  negative,  then  the  penetration  process  will  always  be  complete 
before  the  rod  is  completely  consumed  by  plastic  deformation.  The  final 
undeformed  length  is  obtained  from  the  solution  of 

Uf/L)2a  -  1  -  ap2V2/b(1+u)2  (28) 

The  total  depth  of  penetration  is  given  by  Equation  (26a)  with  given  in 
Equation  (28). 

For  positive  values  of  a,  the  situation  is  somewhat  more  complicated. 

If  b/a  =  2c2/(1+2e)  >  u2V2/(1 +u)2,  then  penetration  ceases  before  the  rod  is 
completely  consumed  by  plastic  deformation.  The  final  undeformed  length  is 
still  given  by  Equation  (28)  with  the  total  depth  of  penetration  again  given 
by  Equation  (26a). 

The  only  remaining  possibility  is  a>0  with  b/a  *  2c2/(1+2e)  < 
ii2V2/(l  +p)2.  In  this  case,  the  rod  will  be  completely  consumed  by  plastic 
deformation  before  penetration  ceases.  The  total  penetration  depth  is  then 

zf  -  L/u  (29) 

according  to  Equation  (26a). 

The  other  case  in  which  Equation  (8)  is  linear  in  u  is  for  y2  -  1.  In 
this  case 

u  -  (1/2)  v  +  (1 -x2)c2/v 

For  u  to  be  positive  requires  a  rod  velocity  such  that 


(30) 


Condition  Equation  (31)  will  always  be  satisfied  for  x 2  <  1  until  v 
decreases  to  zero.  For  A2  >  1  there  is  a  critical  initial  striking  speed 
[2(A2~1)]1/2  c.  For  actual  initial  impact  velocities  less  than  this  value, 
penetration  will  not  occur.  For  impact  speeds  greater  than  this  critical 
value,  penetration-with-consumption  will  occur  and  will  continue  until  the 
rod  is  decelerated  to  [2(X2~1 )]1  ^2  c  when  penetration  will  cease. 

For  u  to  be  smaller  than  v  in  this  case  (u2  *  1 )  requires  that 

(1/2)v2  -  (1-X2)c2  >  0  (32) 

Here  the  inequality  will  always  be  satisifed  for  X2  >  1  but,  from  the 
discussion  above,  this  will  include  situations  of  zero  penetration.  For  X2 
<  1  there  is  a  critical  initial  striking  speed  [2(1 -A2)]1 /2  c  below  which 
rigid  body  penetration  occurs  and  Equation  (17)  applies.  For  initial 
striking  speeds  above  this  critical  value  there  will  be  some  penetration- 
with-consumption  until  the  rod  is  decelerated  to  [2(1-A2)J1/2  c  after  which 
some  further  rigid  body  penetration  occurs,  for  which  Equation  (16)  applies. 

When  u2  W  1  and  X2  W  1  Equation  (8)  is  quadratic  in  u  and  has  the 
solution 

u  -  (v  -  Cu2v2  -  2(1-u2)(1-x2)c2]1/2}/(1-u2)  (33) 

Whenever  v  >  [2(1  -p2)(1  -X2)]1  /2c/y,  the  root  in  Equation  (33)  is  real  and 
the  rod  penetrates  the  target.  The  complexity  of  this  relationship  between 
u  and  v  precludes  the  possibility  of  finding  any  exact  solutions  to  the 
differential  equations  beyond  that  given  in  Equation  (7).  The  system  is 
thus  integrated  numerically. 

5.  LIMIT  AT  HxGH  IMPACT  SPEEDS 

For  v2,  sufficiently  large  in  comparison  to  c2,  Equation  (33)  gives  the 
approximate  result 

u  =  v/(  1  +p)  (3*0 
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which  is  identical  to  Equation  (18).  In  this  case  Equation  (19)  applies, 
and  Equations  (26)  and  (26a)  follow.  Rearranging  Equation  (26a)  gives 


uzf/L  -  1  -  lf/L  (26b) 

At  the  highest  impact  velocities  the  penetrator  is  completely  consumed 
(if  »  0)  and  from  Equation  (26b)  it  follows  that  the  limit  of  yZf/L  is  unity 
as  V/c  becomes  very  large.  This  is  obvious  from  the  numerical  calculations 
shown  in  the  figures. 
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SECTION  III 


RESULTS  AND  DISCUSSION 

Figure  3  shows  some  calculated  results  for  cases  favorable  to  penetra¬ 
tion.  In  these  cases,  by  comparison  to  the  target  the  penetrator  is 
stronger  (A2  »  1/3)  and  more  dense  (y2  =  0.6-4 ).  This  figure  shows  the 
effect  of  one  element  introduced  by  the  present  theory,  viz.,  the  mushroom 
strain  at  the  projectile  end.  As  can  be  seen  from  Equations  (9)  and  (9a) 
this  strain,  which  is  negative,  directly  increases  the  projectile  decelera¬ 
tion.  The  obvious  physica'l  interpretation  is  that  the  incoming  projectile 
momentum  is  resisted  by  a  larger  target  area.  For  zero  strain  the  results 
of  the  present  analysis  are  nearly  identical  to  those  calculated  from  Tate’s 
analysis.  However,  as  the  impact  end  strain  is  assigned  increasingly  larger 
magnitudes,  the  penetration  depth  is  greatly  reduced.  It  must  be  concluded 
that  this  mushroom  strain  plays  a  highly  significant  role  in  the  penetration 
process. 

Figure  3  shows  that  for  strains  smaller  in  magnitude  than  0.5  the 
penetration  at  first  increases  with  impact  velocity,  reaches  a  maximum,  and 
then  decreases.  This  is  observed  in  the  calculations  based  upon  Tate's 
theory  also,  and  has  been  discussed  at  length  in  his  1969  paper  (Reference 
2).  Based  upon  the  present  analysis,  it  seems  likely  that  when  the  proper 
dependence  oi  e  upon  V  is  discovered,  penetration  depth  will  be  predicted  to 
increase  monotonlcaily  with  impact  speed. 

Figure  U  shows  some  calculated  results  for  cases  in  which  the  penetra¬ 
tor  and  target  are  of  identical  materials.  They  have  equal  strengths  (A2  « 
1)  and  densities  (y1-  »  1).  This  figure  shows  the  effect  of  the  second 
major  element  introduced  by  the  present  theory,  viz.,  the  relative  velocity 
term  in  the  equation  of  motion.  This  is  most  easily  seen  in  Equation  (9a). 


Normalized  Impact  Velocity  for 


Here  the  term  ~cd  that  occurs  in  Tate's  equation  is  reduced  by  the  term 
(1/2  +  e)(v  -  u)2  in  the  present  theory.  For  e  =  0,  Figure  4  shows  that 
substantially  larger  penetrations  are  calculated  from  the  present  theory 
than  from  the  Tate  theory  over  the  range  of  impact  velocities  bridging 
between  rigid  body  penetration  and  the  limiting  case  at  high  velocities. 

This  effect  is  not  so  strong  as  the  mushroom  strain  effect,  however,  and  the 
figure  shows  a  reverse  to  smaller  penetration  depths  at  two  non-zero 
stra ins. 

Figure  5  shows  some  calculated  results  for  cases  unfavorable  to 
penetration.  In  these  cases,  by  comparison  to  the  target  the  penetrator  is 
weaker  (A2  =  3)  and  less  dense  (p2  =  1.44).  This  figure  shows  that  for  e  = 
0,  calculations  based  upon  the  present  theory  are  virtually  identical  to 
those  based  on  Tate's  theory.  But  there  remains  a  highly  significant  effect 
due  to  impact  face  strain. 


SECTION  IV 


CONCLUSIONS 

The  present  analysis  provides  a  one-dimensional  penetration  theory  that 
is  more  accurate  than  its  antecedent. 

Two  new  effects  are  introduced  in  the  present  theory:  the  change  of 
momentum  of  material  crossing  the  plastic  interface,  and  the  strain  discon¬ 
tinuity  at  this  interface.  Of  these  two,  the  latter  has  the  greater  impact 
on  calculated  results. 

The  present  analysis  is  capable  of  describing  impact  penetration  data 
without  recourse  to  abnormally  large  strength  parameters  for  the  target  and 


penetrator  materials. 
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